Subharmonic gap structure in superconducting scanning tunneling microscope 

junctions 
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We observe a subharmonic gap sturucture (SGS) and the Josephson effect in superconducting 
scanning tunneling microscope junctions with resistances below fOO kfl. The magnitude of the 
n — 2 SGS is shown to scale with the square of the junction normal state conductance, in agreement 
with theory. We show by analyzing the Josephson effect in these junctions that the superconducting 
phase dynamics are strongly affected by thermal fluctuations. We estimate the linewidth of the 
Josephson oscillations due to phase fluctuations, a quantity that may be important in modern 
theories of the subgap structure. While phase fluctuations may smear the SGS current onsets, we 
conclude that the sharpness of these onsets in our data is not limited by fluctuations. 
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The enhancement of the subgap tunneling current in 
superconducting (SC) tunnel junctions at voltages eV — 
2A / n, where A is the SC gap and n is an integer, was first 
observed many years ago [l| . The phenomenon is termed 
the subharmonic gap structure (SGS) and is character- 
ized by temperature independent excess current onsets 
at these voltages. Renewed interest in this effect over 

the last decade HUSH, 0. B B D3. El 03 was 

driven in part by advances in fabrication of controllable 
superconducting point contacts. The emergence of new 
systems, such as Josephson junction qubits where quasi- 
particle currents play a role in the decoherence of the 
qubit state EHj al so calls for a better understanding of 
subgap quasiparticle tunneling processes. 

Early theoretical attempts to describe this effect in 
"good" (low transparency) tunnel junctions offered sev- 
eral possible mechanisms, most notably multiparticle 
tunneling E3 (MPT) and Josephson self coupling EH- 
Both mechanisms can account for the observed subhar- 
monic structure, but their predictions for the dependence 
of the effect on the junction transparency is different Ell- 
in the limit of high transparency contacts, Klapwijk et 
al. E3l proposed multiple Andreev reflections (MAR) as 
a plausible mechanism. Experiments using controllable 
atomic break junctions 0,13 have traced the behavior of 
the SGS through a wide range of junction transparencies 
from the tunneling to the point contact regimes. These 
have led to a theoretical effort 0, 0, to interpolate 
between the descriptions of SGS in the low (MPT) and 
high (MAR) transparency limits. 

The theories of Refs. 0, Q, while successful in ex- 
plaining the experimental results, formally rely on a well 
defined time evolution of the SC phase across the junction 
with the Josephson frequency to — 2eV/ h in the presence 
of a bias voltage. Indeed it is well known that the SGS 
in atomic junctions is smeared in the presence of rf noise 
El- However, a theoretical description of this smearing 
within the framework of Refs. and its connection 

with noise induced phase fluctuations and broadening of 
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FIG. 1: Typical tunneling conductance for a superconducting 
STM vacuum tunnel junction at T=2.1 K and high junction 
resistance Rm ~ 10 Mf2. 



the Josephson oscillations is not yet established. 

We report here on experiments using SC scanning tun- 
neling microscope (STM) tips to form vacuum tunnel 
junctions against SC films. For moderate junction trans- 
parencies [lfl we observe the Josephson effect and the 
subgap structures at eV — 2 A/2 and eV = 2 A/3. In 
contrast with the break junction experiments |?l lil floj| . 
where thermal and noise induced phase fluctuations were 
kept minimal and their effect on the SGS was neg- 
ligible, the noise temperature in the present work is 
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quite high. With the observation of a smeared SGS in 
the current-voltage characteristics of our junctions, and 
a measurement of the noise temperature based on the 
Josephson effect, this work represents a first attempt at 
studying the effects of phase fluctuations on the SGS, 
which so far has been only qualitative. Furthermore, this 
work emphasizes the onset of the SGS in the low trans- 
parency tunneling regime, where most of the existing lit- 
erature is concerned with high transparency near-contact 
junctions. 

All measurements were performed at T=2.1 K in a 
home built STM. The scan signals, thermometer, and 
bias voltage lines were filtered before entering the cryo- 
stat using commercial feedthrough filters. Additional fil- 
tering of the bias leads was provided by an RC filter 
mounted near the junction at 2.1 K. The tunneling cur- 
rent is measured with a current-voltage converter (1VC) 
at room temperature that is connected to the junction 
via a coaxial cable. A resistor mounted at the input of 
the IVC and the capacitance of the cable act as a low 
pass filter above ~20 MHz. Since we do not employ 
special microwave filters in these experiments, and the 
effectiveness of the filters that we use severely degrades 
at microwave frequencies, we expect high frequency noise 
from the room temperature electronics to reach the junc- 
tion and contribute to an elevated noise temperature. 

Superconducting STM tips were prepared as detailed 
in Ref. H2 by depositing a Pb(5000 A)/Ag(30 A) prox- 
imity bilayer onto mechanically cut Pto.slro.2 tips. A 
freshly cleaved graphite crystal was placed near the tips 
at the time of the deposition and a Pb/Ag film was de- 
posited on it as well, serving as the sample. We verified 
the quality of the tip by scanning the tip over the sample 
in the imaging mode. Typical spatial resolution obtained 
with these tips is ~10 A. Current-voltage (J— V) and con- 
ductance (dl/dV) curves were then taken to verify the 
superconducting properties of the junctions. A typical 
dl/dV curve is shown in Fig. ^ where the sharp conduc- 
tance peaks at voltages V — ±2Apb/e, the small leakage 
conductance below the gap, and the Pb phonon structure 
confirm that these are good tunnel junctions. 

A sequence oi I — V curves was recorded as the junc- 
tion normal state resistance was lowered by moving the 
tip closer to the sample, thereby increasing the junc- 
tion's transparency as |T| 2 cx Rff . Fig.|2]shows selected 
I—V curves from such a sequence, for junction resistances 
ranging from 70 kSl to 8 kfi. New features in the tun- 
neling characteristics emerge as the junction resistance is 
lowered - a current peak near zero bias that is the sig- 
nature of pair tunneling |22j . and the subharmonic gap 
structure at voltages 2A/2e and 2A/3e. When the data 
are numerically differentiated (Fig. inset) the subgap 
structure is seen as conductance peaks at V=f.35 mV 
and V=0.85 mV for the n = 2 and n = 3 subharmonics 
respectively. 

In order to determine the transparency dependence of 
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FIG. 2: I— V characteristics for Rn between 70 kfi (squares) 
and 8 kQ (diamonds) showing subharmonic gap structure and 
fluctuation dominated Josephson effect. The magnified (J x 5) 
portion of the lowest resistance curve shows the n — 3 struc- 
ture more clearly. Inset: numerically smoothed and differen- 
tiated data. 

the SGS magnitude, we plot the peak conductance, af- 
ter subtracting the single quasiparticle background, for 
the n — 2 structure as a function of the junction normal 
state conductance on a logarithmic scale in Fig. |21 The 
linear relation obtained from this graph suggests that 
Gsgs = Gp ca k - Gsingic oc G% , and from a fit we ob- 
tain a = 2.3 ± 0.2, in good agreement with the expected 
|T| 4 ~ G 2 N dependence. We observe a similar dependence 
in all junctions measured to date, with the exponent a 
ranging from 1.8 to 2.3. Strictly speaking, it is the cur- 
rent jump at this voltage that is expected to scale with 
|T| 4 , but since the width of this onset in our data is inde- 
pendent of Gat, the behavior of the differential conduc- 
tance peak reflects that of the current jump. Unfortu- 
nately, the present data do not allow us to determine the 
transparency dependence of the n = 3 feature, although 
its magnitude relative to the n = 2 feature is consistent 
with the theoretically predicted |T| 6 ~ G Z N dependence. 

Were the junctions single channel contacts, their trans- 
parency would have saturated to I at Rm — 12.9 kfi. The 
fact that we still observe tunneling-like / — Vs down to 
8 kil suggest that our tips accomodate a few channels. 
This conclusion is also consistent with the ~10 A spatial 
resolution that we estimate from working in the imaging 
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FIG. 3: The differential conductance at eV = 2A/2 versus the 
junction normal state conductance (log scale) for two different 
junctions (sample 1 is shown in Fig.|2Jl. The single quasipar- 
ticle background that we have subtracted was deduced from 
curves measured at higher junction resistances for which the 
normalized conductance is independent of Rn- The solid line 
represents a slope of 2.3. 

mode. The relatively large scatter in the data in Fig. [3] 
can thus be understood as resulting from atomic rear- 
rangements in the tip during the time of the experiment. 
Such rearrangements may change the transparencies of 
the individual channels in the junction; while there are 
many different combinations of channel transparencies 
that can result in the same normal state conductance, 
the I — V characteristics below the SC gap are highly 
nonlinear and therefore sensitive to the particular chan- 
nel content of the junction 0- 

In light of the theoretical descriptions of the SGS 
m Refs. 0, H @ it is important to further character- 
ize the phase dynamics of these junctions. The evolu- 
tion of the Josephson current as a function of the junc- 
tion normal state resistance allows us to quantify the 
strength of the phase fluctuations in the junction. Be- 
cause of the junctions' small size and the high resis- 
tance associated with it, the Josephson binding energy 
Ej = 7rfi.A/4e 2 i?Ar, which is the energy scale for the 
coupling of the phases across the junction, is compa- 
rable to ksT. In this regime, where Ej m knT, the 
motion of the phase difference is diffusive [23I |24| , and 



FIG. 4: A plot of \J (4e/h)A/V p vs. the junction normal state 
conductance. Inset: fits of the I — V characteristics using the 
phase diffusion model of Ref. |23| . 



the pair current is therefore associated with finite volt- 
ages. The I —V characteristics due to pair tunneling for 
small Ej /UbT were shown by Ivanchenko and Zil'berman 
H3 to have the form I(V) = A x V/(V 2 + V p 2 ), where 
A = I 2 c Z cm /2, I c = 2eEj/H, V p = {2e/h)Z env k B T n , Z cav 
is the impedance of the junction's environment, and T n 
is the effective noise temperature. We fit the I— V curves 
using the equation above with A and V p the only fitting 
parameters (Fig. 0] inset). A plot of yj (Ae/h)A/V p as 
a function of Gat yields a linear graph with a slope of 
I c Rn I W^bTu [23|. Such a plot (Fig. is constructed 
from analysis of the data of Fig. |2J From the slope of the 
graph, and using the Ambegaokar-Baratoff relation 
to find I c Rn — 1-67 mV in this junction, we estimate the 
noise temperature T n — 33 ± 1 K |26| , and the environ- 
ment impedance Z cnv = 78 ± 7f2 close to the impedance 
of free space as expected [23 ■ 

The quantity of importance to the dicussion here is the 
linewidth of Josephson oscillations around <j> = 2eV/h. 
The mechanism for the formation of the SGS within the 
theory of Bratus' et al. 3 is the inelastic Andreev scat- 
tering of quasiparticles impinging on the junction into 
n sidebands. The splitting in the spectrum of scat- 
tered waves into sidebands is induced by the time vary- 
ing phase across the junction with the Josephson fre- 
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FIG. 5: Compaison of the widths of the n = 2 SGS (top 
panel) and the gap edge (bottom panel). Note the different 
conductance scale, and the offset voltage scales. 



quency w — 2eV/h. Our experiments show that in the 
presence of thermal phase fluctuations there is a con- 
siderable broadening in the time evolution of the phase 
[22| . where the linewidth T around the Josephson fre- 
quency depends on the effective noise temperature and 
the impedance of the junction's environment. In the 
RSJ model this broadening can be estimated [2^] as 
T w (2e/fi) 2 Z onv fcsT„ = 2eV p /h. For the data shown 
here V p = 0.11 ± 0.01 mV, resulting in T ~ 300 GHz. 

We expect that the effects of phase fluctuations would 
be to broaden the subgap features 0,0] similarly to the 
broadening of microwave induced Shapiro steps |23 • It is 
therefore important to determine whether the broadening 
of the SGS features in our junctions (Fig. |2 inset), with 
full width of approximately 0.35 mV, arise from phase 
fluctuations. Since any subgap quasiparticle tunneling 
process would eventually involve the injection of a quasi- 
particle at the gap edge of the SC electrode, the SGS 
current onsets can be no sharper than the current rise 
at eV — 2A. A well documented phenomenon in Pb, of 
which our junctions are made, is the apparent broaden- 
ing of the gap edge in tunneling experiments 

minim. 

This broadening is intrinsic in bulk Pb, and is a result of 
the ~10% anisotropy in the Pb order parameter [32ll33j |. 
Figure shows a comparison between the conductance 
peaks at the gap edge and at the n — 2 SGS. The 
fact that the widths of these peaks are approximately 



the same, and larger than the characteristic broadening 
of the Josephson oscillations (h/2e)T, suggests that the 
sharpness of the SGS in our junctions is not fluctuation 
limited, but rather it is intrinsic to the superconductors 
comprising the junction. We would expect that with in- 
creasing strength of the phase fluctuations, the width of 
the SGS would eventually increase above the width of 
the gap edge. 

To summarize, we have observed the n = 2 and n = 3 
subharmonic gap structure in superconducting STM vac- 
uum junctions, and found a G 2 N dependence for the mag- 
nitude of the n — 2 SGS. This experiment has been 
done under conditions where strong phase fluctuations 
are present, and by analyzing the Josephson effect in 
these junctions we have determined the broadening of the 
Josephson oscillation linewidth due to these fluctuations. 
Since modern theories for the origins of the excess sub- 
gap current 0,H,IE0| suggest that the subharmonic gap 
structure and the ac Josephson effect are intimately re- 
lated, it could be expected that strong phase fluctuations 
will alter the shape of the SGS. Nevertheless, we conclude 
that the apparent smearing of the observed SGS in our 
data is most likely due to the intrinsic gap anisotropy 
in Pb, rather than phase fluctuations. Experiments such 
as described here, could in principle, provide a testing 
ground for theoretical descriptions of the effects of fluc- 
tuations on the SGS. The complications introduced by 
the intrinsic width of the gap edge in Pb, can be easily 
overcome by fabricating the junctions from Pbi-^Bi^ al- 
loys with x ~ 0.10 — 0.15 in which the mean free path 
is shorter, and as a result, the gap edge is sharper, with 
characteristic widths at 2 K of 0.068 meV and 0.048 meV, 
for x = 0.10 and 0.15 respectively HlUll. 

This experiment demonstrates that STM studies with 
well characterized SC tips are suitable for careful mea- 
surements of the subgap structure. Additional advan- 
tage of STM experiments is that in contrast with break 
junctions that are necessarily symmetric, the STM setup 
offers the possibility to study asymmetric junctions (e.g. 
Pb/I/NbSe^) [13. We note that the STM has been used 
before 0, to observe the SGS, but unfortunately a 
quantitative analysis of the phase fluctuations in those 
experiments was not reported. We argue that the mea- 
surement of the phase fluctuations rather than their sup- 
pression may further help to test the existing theories. 
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